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PREFACE 


This  Memorandum  describes  a  computer  program  that  demonstrates 
the  effect  of  denying  the  use  of  a  portion  of  a  highway  transporta¬ 
tion  network.  This  application  utilizes  previous  RAND  research  on 

* 

network  flows  and  on  highway  capacity.  The  program  is  being  used  in 
the  continuing  RAND  research  on  the  effective  employment  of  air  power, 
including  the  employment  of  tactical  air  forces. 

The  model  was  originally  intended  for  internal  use  at  RAND,  but 
a  number  of  other  agencies  have  indicated  that  it  might  prove  useful 
in  their  research,  including  the  Weapons  Systems  Evaluation  Croup  (WSEG) 
and  the  Operations  Analysis  Office,  Hq  USAF  (AFGOA) .  It  should  be  of 
interest  both  to  those  concerned  with  targeting  strikes  against  road 
networks  and  to  those  concerned  with  allocating  road  repair  and  improve¬ 
ment  efforts. 


L.  R.  Ford,  Jr.,  and  D.  R.  Fulkerson,  Flows  in  Networks,  The 
RAND  Corporation,  R-375-PR  (DDC  No.  AD  287499),  August  1962  (also 
published  by  Princeton  University  Press,  1962);  R.  D.  Wollmer,  Some 
Methods  for  Determining  the  Most  Vital  Link  in  a  Railway  Network,  The 
RAND  Corporation,  RM-3321-ISA,  April  1963;  L.  P.  Holliday,  A  Method 
for  Estimating  Road  Capacity  and  Truck  Requirements  (U)  ,  The  RAND 
Corporation,  RM-3331-ARPA,  November  1963  (Confidential). 
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SUMMARY 


This  Memorandum  describes  a  computer  program  designed  to  evaluate 
the  capability  of  a  transportation  network  to  deliver  supplies  to 
destinations  as  road  segments  or  arcs  making  up  the  network  are 
successively  destroyed  and  repaired.  The  program,  written  in  FORTRAN  IV, 
can  be  adapted  easily  for  use  on  any  of  several  large-scale  computers. 

As  inputs,  the  program  requires  a  description  of  the  transportation 
road  system  under  consideration,  and  of  the  cargo-carrying  vehicles 
operating  on  the  network.  Given  the  basic  data,  the  program  furnishes 
a  profile  of  maximum  cargo  flow  as  a  function  of  the  number  of  vehicles 
made  available  to  the  system,  and  then  selects  and  destroys  that 
vulnerable  link  in  the  network  which  reduces  the  cargo  flow  rate  most 
severely.  The  program  repeats  these  steps  until  flow  on  the  network 
is  totally  stopped  or  the  predesignated  number  of  links  have  been 
destroyed.  The  program  then  steps  to  the  next  "day"  or  "period," 
restores  to  service  all  previously  destroyed  links  that  have  been 
repaired  by  this  date,  and  repeats  the  process  of  profile  generation 
and  link  removal. 

* 

This  model  uses  the  Fulkerson  Out-of-Kilter  Algorithm 
to  generate  the  cargo  flow  profile  as  a  function  of  vehicles  in  the 

•kk 

system,  and  uses  an  algorithm  of  Wollmer  to  determine  the  single 
most  critical  link  in  a  network. 

The  program  presently  will  accept  a  network  of  up  to  1000  links, 
but  this  number  may  be  modified  to  suit  the  capacity  of  a  particular 
computer  system.  By  properly  describing  the  network,  some  combinations 
of  rail,  road,  and  water  transportation  can  be  analyzed.  The  program 
should  be  a  useful  tool  in  targeting,  in  logistics  system  analysis, 
in  allocating  funds  for  expansion  of  transportation  systems,  and  allocating 
road-repair  efforts. 


Ford  and  Fulkerson,  op  cit.;  D.  R.  Fulkerson,  An  Out-of-Kilter 
Method  for  Minimal  Cost  Flow  Problems,  The  RAND  Corporation,  P-1825, 
April  1960  (also  published  in  J.  Soc,  Ind,  Appl.  Math.,  Vol.  9, 

March  1961,  pp.  18-27). 

** 

Wollmer,  op.  cit. 
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I.  INTRODUCTION 


The  general  tactical  interdiction  problem  is  to  select  target 
elements  from  among  the  vehicles,  road  links,  rail  nets,  waterways, 
and  air  lines  supplying  a  combat  force  in  order  to  effectively  reduce 
the  combat  capability  of  that  force.  In  order  to  select  an  optimal 
set  of  targets  to  attack,  a  general  interdiction  model  would  presumably 
consider  both  the  weapons  and  alternatives  available  to  the  interdictor, 
and  the  current  configuration  and  recuperative  ability  of  the  trans¬ 
portation  system.  This  selection  would  change  over  time  as  lucrative 
concentrations  of  vehicles  were  observed,  as  critical  segments  of  the 
network  were  repaired  and  returned  to  service,  and  as  the  weapons  available 
to  the  interdictor  changed.  The  obvious  information  one  might  require 
from  such  a  model  would  be  the  ability  of  the  transportation  system 
to  support  the  combat  force  at  various  levels  of  attack,  or  the  level 
of  attack  required  to  degrade  combat  performance  to  a  specified  value. 

The  more  restricted  situation  considered  by  the  model  described 

here  explicitly  includes  only  a  single  type  of  cargo-carrying  vehicle 

moving  over  a  highway  network.  The  restriction  to  a  single  type  of 

cargo  vehicle  was  imposed  since  it  allows  the  use  of  very  efficient 

algorithms  to  calculate  maximum  flow  through  a  network,  and  it  was 

decided  to  consider  highway  rather  than  water,  air,  or  rail  transporta- 

* 

tion  since  research  by  L.  P.  Holliday  allows  the  network  concept  of 
arc  capacity  to  be  applied  to  highways.  The  key  assumption  in  this 
model  is  that  the  maximum  speed  with  which  vehicles  can  flow  across 
an  arc  is  independent  of  the  volume  of  flow  on  the  arc. 

Given  this  assumption,  the  highway  transportation  system  is  described 
as  a  network  of  nodes  and  directed  arcs  together  with  the  vehicles 
that  move  on  these  arcs.  Nodes  may  be  towns,  intersections,  or  any 
points  at  which  it  is  convenient  to  distinguish  between  the  road 
characteristics  on  either  side  of  the  node.  Each  arc  in  the  network 

* 

Holliday,  op.  cit. 
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is  characterized  by  the  name  of  the  node  at  which  it  originates,  the 
name  of  the  node  at  which  it  terminates,  the  maximum  flow  capacity  on 
the  arc  in  vehicles  per  time-unit,  the  number  of  time-units  required 
for  a  vehicle  to  move  across  the  arc,  and  the  number  of  time-units 
required  to  repair  the  arc  if  it  should  be  rendered  impassable.  If  an 
arc  is  assigned  a  repair  time  of  zero,  it  is  considered  invulnerable 
and  will  not  be  cut.  Arcs  may  have  minimum  flow  requirements  assigned 
in  tons  per  day,  permitting  demands  at  various  destinations  to  be 
specified. 

Other  inputs  to  the  model  are  the  number  of  vehicles  in  the 

inventory,  vehicle  tonnage  capacity,  vehicle  in-commission  rate,  en 

route  stop  factor  (to  convert  en  route  time  to  total  time),  the  number 

of  arcs  that  can  be  cut  per  day,  the  number  of  days  the  model  is  to 

consider,  and  the  conversion  factors  between  time-units,  hours,  and 

* 

operating  days.  Holliday  discusses  the  methodology  for  determining 
road  capacities,  maximum  road  speeds,  and  time  adjustments  required 
by  one-way  use  of  the  roads,  waiting  time,  and  convoying. 

Although  the  actual  number  of  vehicles  in  the  inventory  is 
specified  to  the  model,  the  problem  solved  is  that  of  determining 
throughput  capability,  measured  in  vehicles  per  unit  time  or,  equivalently, 
tons  per  day,  as  a  function  of  the  number  of  usable  arcs  in  the  network 
as  the  number  of  vehicles  in  the  system  varies  from  zero  up  to  network 
saturation  level.  This  is  done  so  that  when  use  of  an  arc  is  denied, 
the  resulting  profile  of  throughput  as  a  function  of  vehicles  available 
to  the  system  can  be  compared  directly  to  the  previous  profile  to 
determine  the  "vehicle-worth"  of  the  newly  removed  arc. 

This  model  makes  repeated  use  of  the  Fulkerson  Out-of-Kilter 
algorithm,  which  is  an  efficient  algorithm  for  constructing  minimal 
cost  flows  in  networks  in  which  the  arcs  have  costs,  and  both  upper 
and  lower  bounds  on  permissible  flow. 

In  selecting  the  most  critical  arc  for  removal,  the  model  uses 
an  algorithm  developed  by  R.  D.  Wollmer.  The  essence  of  this 
algorithm  is  that  in  considering  a  network  with  a  maximal  flow,  F  , 

★ 

Ibid. 

Ford  and  Fulkerson,  op.  cit.;  Fulkerson,  op.  cit. 

Wollmer,  op.  cit. 
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given  the  usable  arcs  of  the  network  and  the  actual  number  of  vehicles 
available  to  the  system,  there  is  no  need  to  consider  those  arcs  with 
zero  flow  as  candidates  for  removal.  Removing  such  arcs  would  con¬ 
tribute  nothing  toward  the  goal  of  identifying  that  arc  which,  if 
removed,  would  most  drastically  reduce  throughput. 
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II,  DESCRIPTION  OF  THE  MODEL 


The  problem  of  determining  maximum  flow  through  a  network  as  a 
function  of  the  number  of  vehicles  in  the  network  is  seen  clearly 
when  formulated  as  a  parametric  linear  programming  problem. 

Let  t. .  =  time-units  required  for  a  vehicle  to  move  from 

i  j 

Node  i  to  Node  j 

u.  .  =  maximum  flow  capacity  in  vehicles  per  unit  time  on  the 
1J  arc  (i ,  j) 

l  .  =  minimum  flow  requirement  in  vehicles  per  unit  time  on 


i  j 


the  arc  (i,j) 


x. .  =  number  of  vehicles  per  unit  time  passing  over  the  arc 

1J  <i,j) 

d^  =  minimum  number  of  vehicles  per  unit  time  required  at 
Destination  i 

T  =  number  of  vehicles  available  to  the  system. 

Connect  an  artificial  source,  Node  0,  to  all  true  sources,  and 
an  artificial  destination,  Node  n+1,  to  all  true  destinations.  Also 
introduce  an  arc  directed  from  the  artificial  destination  to  the 
artificial  source.  Arcs  connected  to  artificial  nodes  are  generally 
assigned  no  cost  and  infinite  capacity.  The  problem  is  then 


(1) 


maximize  x 
x.  . 


n+1 ,0 


subject  to 


(2) 


£  x,  .  -  £  x..  =  0  for  each  j 
i  i  Ji 

*  xtj  ^  for  all  arcs  (i,j) 


where 


d^  for  all  destinations  i,  and 
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(3) 


£  t. .x  . 
i  iJ  ij 


<  T. 


While  the  problem  is  seen  clearly  when  formulated  in  this  way,  the 
presence  of  Constraint  (3)  precludes  solution  by  standard  capacitated 
transportation  algorithms,  while  the  sice  of  realistic  networks  pre¬ 
cludes  solution  by  standard  linear  programming  algorithms.  The  model 
therefore  considers  the  inverse  problem  --  that  of  plotting  the  number 
of  vehicles  required  by  the  system  as  a  function  of  the  number  of 
vehicles  flowing  through  the  system.  In  the  previous  notation  this  it. 

(4)  minimize  £  t. .x. .  , 

x  iJ 

ij 


with 


(5) 


£  x. .  -  £  x..  =0  for  each  j 
i  U  i  J1 


\  j l.  .  <  x.  .  <  u.  .  for  all  arcs  (i,  i)  , 
ij  ij  ij 


where  flow  in  the  system  is  successively  incremented  by  incrementing 
the  lower  bound  on  the  return  flow,  £n+-^  q  •  Solving  a  sequence  of 
such  problems  yields  points  on  the  profile  of  vehicles  required  by 
the  system  as  a  function  of  vehicles  flowing  through  the  system.  This 
is  termed  a  flow  profile.  The  number  of  such  profile  points  to  be 
generated  is  specified  as  an  input  to  the  program. 

Each  time  a  flow  profile  is  to  be  generated,  the  model  first 
computes  the  maximum  flow  through  the  network  constrained  only  by  the 
arc  capacities  and  the  demands  at  the  various  destinations.  This 
maximum  flow  is  divided  by  the  number  of  requested  profile  points  to 
obtain  the  quantity  of  flow  by  which  the  system  will  be  successively 
incremented.  The  vehicle  minimization  problem  stated  in  Expressions 
(4)  and  (5)  is  actually  solved  at  each  level  of  flow,  F,  and  also 
at  F+l,  which  approximates  determining  the  derivative  of  vehicle 
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requirements  at  F,  and  facilitates  plotting  the  profile.  If  the 
actual  number  of  vehicles  declared  to  be  available  to  the  system  falls 
between  the  values  computed  at  two  successive  profile  points,  a  series 
of  interpolations  are  carried  out  that  yield  a  profile  point  within 
50  vehicles  of  the  actual  number.  If  the  flow  that  the  model  is  attempt 
ing  to  force  through  the  system  is  less  than  the  sum  of  all  the  demands 
imposed  on  the  destinations,  the  model  reduces  the  demands  to  the 
proper  proportion  of  the  original  demands.  If  the  model  determines 
that  either  the  actual  or  the  proportionately  reduced  demands  cannot 
be  satisfied,  no  further  profiles  are  generated  for  the  current  day. 

The  network  is  assumed  to  be  interdicted. 

After  the  profile  is  generated  for  a  fixed  network  configuration, 
the  program  selects  an  arc  to  remove.  First,  a  maximum  flow  capacity 
is  imposed  on  the  return  arc  (n+1,0),  which  prevents  maximum  flow  in 
the  network  from  exceeding  that  flow  which  can  be  sustained  by  the 
actual  number  of  vehicles  declared  available  to  the  system.  Vulnerable 
arcs  are  then  removed  in  turn  while  an  attempt  is  made  to  maximize  flow 
through  the  network  subject  to  the  new  constraint  on  maximum  flow.  As 
flow  is  maximized  with  successive  arcs  removed  for  test,  all  arcs  on 
which  there  is  zero  flow  are  removed  from  further  consideration.  That 
arc  is  finally  selected  for  removal  which  causes  maximum  flow  through 
the  network  --  constrained  by  the  actual  number  of  vehicles  available  -- 
to  be  minimized.  If  several  arcs  yield  the  same  minimum  flow  value, 
f,  that  arc  is  removed  which  requires  the  greatest  number  of  vehicles 
to  be  used  at  the  flow,  f.  Arc  removal  is  accomplished  in  the  program 
by  setting  arc  maximum  capacity  to  zero.  At  the  time  the  arc  is 
removed,  its  return-to-service  time  is  updated  by  the  arc  repair  time. 
When  the  return-to-service  time  is  reached,  the  arc  is  restored  to 
full  use. 


III.  PROGRAM  SEQUENCE  OF  OPERATION 


The  program  operates  in  daily  cycles  for  the  number  of  days  spec¬ 
ified  by  the  user.  At  the  start  of  Day  k,  all  arcs  with  appropriate 
return-to-service  dates  are  restored  to  full  capacity.  An  attempt  is 
then  made  to  generate  the  flow  profile  at  the  start  of  Day  k.  If  an 
insufficient  number  of  arcs  have  been  repaired  by  this  date,  and  no 
flow  is  possible,  the  message  "FL0W  T0 TALLY  ST0PPED  0N  DAY  K"  will  be 
printed  and  the  program  will  cycle  to  Day  k-!-l.  If  flow  through  the 
network  is  possible,  but  not  in  a  pattern  appropriate  to  meet  demands 
imposed  at  various  destinations,  the  message  "MIN  DEMANDS  AT  S0ME  DES- 
TINATI0N  CANN0T  BE  SATISFIED"  will  be  printed,  and  the  program  will 
cycle  to  Day  k+1.  If  neither  of  these  conditions  occurs,  the  flow 
profile  at  the  start  of  Day  k  will  be  printed. 

Next,  an  arc  will  be  selected  for  removal,  its  capacity  will  be 
set  equal  to  zero,  and  the  flow  profile  resulting  from  the  new  network 
configuration  will  be  generated  and  printed.  The  cycle  of  arc  removal 
and  flow  profile  generation  continues  until  either  minimum  demand  can¬ 
not  be  met  at  some  destination,  a  preassigned  number  of  arcs  have  been 
removed,  or  flow  is  totally  stopped.  When  any  of  these  conditions 
occur,  the  model  proceeds  to  Day  k+1  and  the  sequence  begins  again. 

The  program  prints  several  other  self-explanatory  data-editing  messages 
and  statements  describing  the  arcs  that  have  been  cut. 
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IV.  INPUT  DESCRIPTION  AND  PREPARATION 

MODEL  CONTROL  DATA 

1.  Number  of  Vehicles  in  the  Inventory  (NTRKS) .  The  model  uses 
this  information  only  to  place  a  marker  (***)  at  the  appropriate  point 
in  the  flow  profile.  The  flow  profile  is  continued  until  network 
saturation. 

2.  Length  of  System  Evaluation  (NDAYS) .  This  is  the  number  of 
days  for  which  the  model  is  to  continue  network  evaluation.  Since 
the  present  formulation  contains  no  random  elements,  it  is  useless  to 
continue  evaluation  beyond  the  point  at  which  network  behavior  begins 
to  repeat  itself. 

3.  Number  of  Arcs  to  be  Removed  Each  Day  (NARCS) .  This  is  the 
maximum  number  of  arcs  that  will  be  removed  each  day.  If  the  network 
is  interdicted  prior  to  this  maximum  number  of  cuts,  arc  removal  will 
cease. 

4.  Number  of  Profile  Points  (NPPTS) .  The  profile  of  flow  versus 
vehicles  will  be  drawn  through  the  plotted  points.  This  number  of 
points  must  be  specified. 

5.  Name  of  Artificial  Source  and  Destination  (SOURCE. SINK) . 

Two  artificial  nodes  are  required.  One  will  be  the  artificial  source 
and  must  be  connected  to  all  true  sources.  The  other  is  the  artificial 
destination  and  all  true  destinations  must  be  connected  to  it.  An 

arc  must  also  connect  the  artificial  destination  to  the  artificial 
source . 

VEHICLE  AND  OPERATING  DATA 

1.  Tons  per  Vehicle  (TNSTRK).  All  vehicles  carry  the  same  load. 

The  model  actually  evaluates  vehicle  flow  per  time-unit  and  then  converts 
this  to  flow  in  tons  per  day. 

2.  Time-Units  per  Hour  (TUPRHR) .  Arc  data  can  be  expressed  in 
arbitrary  time-units.  TUPRHR  is  the  factor  that  defines  the  number 
of  such  arbitrary  time-units  per  hour. 
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3.  Operating  Hours  per  Day  (HRPRDA) .  Vehicles  operate  only  for 
the  number  of  hours  specified  by  this  input. 

4.  In-Commission  Factor  (FINCOM).  The  in-commission  factor  is 
the  fraction  of  the  total  inventory  available  for  dispatch. 

5.  Stop  Factor  (FSTOP) .  This  is  a  multiplicative  factor  that 
can  be  used  to  increase  the  total  number  of  vehicles  in  the  system  to 
account  for  vehicle  stops,  refueling,  etc.  If  K  vehicles  are  actually 
required  to  provide  a  flow  of  F  vehicles  per  unit  time  at  the  final 
destination,  the  model  reports  that  [K* FSTOP /FINCOM]  vehicles  are 
required  to  deliver  F*TNSTRK*TUPRHR* HRPRDA  tons  per  day  to  the  destination. 

ARC  DATA 

Arcs  are  directed  from  Node  i  to  Node  j.  The  names  of  the  nodes 
at  the  beginning  and  end  of  each  arc  are  required.  If  flow  is  possible 
in  both  directions  between  i  and  j,  an  arc  must  be  entered  from  Node 
j  to  Node  i.  When  an  arc  is  cut  between  i  and  j,  the  reverse  arc  from 
j  to  i  is  also  cut.  If  there  are  multiple  arcs  in  both  directions 
between  i  and  j,  the  appropriate  reverse  arc  is  found  by  searching  for 
the  reverse  arc  with  the  identical  repair  time. 

The  program  will  accept  up  to  1000  arcs  and  up  to  500  nodes. 

The  cost  of  obtaining  a  flow  of  one  vehicle  per  time-unit  from  i 
to  j  is  required.  In  the  ground-transportation  problem  this  cost  has 
been  the  number  of  time-units  required  for  one  vehicle  to  traverse  the 
arc  (i,j)  in  both  directions,  given  that  the  vehicle  is  moving  at  the 
most  efficient  speed  for  the  arc.  The  cost  should  include  waiting  time 
required  by  single-lane  operation  on  narrow  roads. 

Maximum  capacity  for  each  arc  is  required  in  terms  of  vehicles  per 
time-unit  passing  an  average  point  on  the  arc. 

Minimum  demand  at  each  destination  can  be  specified  by  placing  a 
minimum  required  flow  in  tons  per  day  on  the  arc  connecting  the  des¬ 
tination  and  the  artificial  destination. 

Repair  times  of  the  arcs  are  entered  in  time-units.  If  a  repair 
time  is  zero,  the  arc  will  never  be  selected  for  removal. 
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The  artificial  arcs  connecting  artificial  to  true  sources,  true 
to  artificial  destinations,  and  artificial  destinations  to  artificial 
sources  should  have  very  large  upper  capacities,  zero  repair  times, 
and  minimum  flow  requirements  where  appropriate. 


INPUT  DATA  FORMAT 


I  denotes  an  integer  right- justified  in  its  field. 

D  denotes  a  decimal  number  anywhere  in  the  field.  If  no  decimal 
point  is  explicitly  entered,  the  decimal  point  is  understood 
to  be  to  the  right  of  the  entire  field. 

A  denotes  any  combination  of  alphabetic  and  numeric  data. 


First  Card 


Col 

1-10 

Number  of  Vehicles  (NTRKS) 

(I) 

Col 

11-20 

Number  of  Days  (NDAYS) 

(I) 

Col 

21-30 

Number  of  Links  (NARCS) 

(I) 

Col 

31-40 

Number  of  Profile  Points  (NPPTS) 

(I) 

Col 

43-48 

Name  of  Artificial  Source 

(A)* 

Col 

51-56 

Name  of  Artificial  Destination 

(A)* 

Second  Card 

Col 

1-5 

Tons  per  Vehicle  (TNSTRK) 

(D> 

Col 

6-10 

In-Commission  Factor  (FINC0M) 

(D) 

Col 

11-15 

Stop  Factor  (FST0P) 

<D) 

Third  Card 

Col 

1-5 

Operation  Hours  per  Day  (HRPRDA) 

(D) 

Col 

6-10 

Time-Units  per  Hour  (TUPRHR) 

(D) 

★ 

Node  names  are  read  as  six  characters,  including  blanks.  There¬ 
fore  node  names  of  less  than  six  characters  must  be  identically  placed 
in  their  field  whenever  written. 


1 
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Fourth  Card 

Col  1-5  The  word  "READY" 

Fifth  Card 

Any  arbitrary  title  information  in  Cols  1-72.  This 
will  be  printed  on  the  outpuL. 

Sixth  Card 

Col  1-4  The  word  "ARCS" 

Arc  Data  Cards 

Each  arc  data  card  has  the  following  format. 


Col 

7-12 

Name  of  First  Node,  i 

(A) 

Col 

13-18 

Name  of  Second  Node,  j 

(A) 

Col 

21-30 

Cost  of  One  Unit  of  Flow,  t. . 

l 1 

(I) 

Col 

31-40 

-J 

Maximum  Capacity  of  Arc,  u 

(I) 

Col 

41-50 

Minimum  Flow  on  Arc,  X,  . 

(I) 

Col 

61-63 

R' ,air  Time  of  Arc,  r. . 

ij 

(I) 

Next  Card 

Col  1-3  The  word  "END" 

Next  Card 

Col  1-7  The  word  "COMPUTE" 


Figure  1  illustrates  a  hypothetical  network  on  which  arc  capacities 
and  costs  have  been  entered.  Node  A  is  the  source,  the  double-lined 
arcs  are  invulnerable,  and  the  demand  at  Node  F  is  440  tons  per  day. 
Artificial  nodes  and  arcs  are  first  introduced  as  shown,  and  it  is 
assumed  that  each  vehicle  in  the  system  carries  3  tons,  that  the  in¬ 
commission  factor  is  0.75,  that  the  stop  factor  is  1.25,  that  vehicles 
operate  15  hours  per  day,  and  that  all  arc  data  are  expressed  in  terms 
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c  •  cost  (time  units  to  transition  arc) 
u  ■  upper  capacity  (vehicles  per  time  unit) 

Double  line  Indicates  invulnerable  arc. 


Fig.  1  --  A  Hypothetical  Network 
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Fig.  2  —  Input  Data  Describing  Example  Network 
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of  30-minute  time-units.  A  three-day  evaluation  is  desired,  with  two 
arcs  to  be  removed  per  day,  and  three  points  generated  on  each  flow 
profile.  Figure  2  indicates  how  these  data  would  be  entered.  The 
Appendixes  depict  the  interdiction  program  and  its  output. 
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Appendix  A 

INTERDICTION  PROGRAM  OUTPUT 


arbitrary  descriptive  oata  for  example 


ARCS 


ARCS 

FROM 

TO 

COST 

UPPER 

LOWER 

TIME 
TO  REI 

ARTIFS 

A 

0 

999 

-0 

0 

A 

B 

A 

30 

-0 

-0 

A 

C 

15 

8 

-0 

1 

» 

D 

11 

15 

-0 

l 

e 

E 

20 

30 

-0 

65 

c 

E 

5 

8 

-0 

-0 

0 

F 

8 

15 

-0 

95 

E 

F 

1 1 

10 

-0 

-0 

F 

ARTIFO 

-0 

999 

4 

-0 

ART 1FD  ARTIFS 

ENO 

COMPUTE 

THIS  PROBLEM  HAS 

-0 

999 

10  ARCS  ANO 

-0 

-0 

8  NODES 

FLOW 

DAY  l 

PROFILE  PRIOR  TO  FIRST  CUT 

TRUCKS 

TONS/OAY 

TRUCKS 

IN  INVENTORY 

THROUGHPUT 

ACTIVE 

306 

720 

184 

345 

810 

207 

460* •• 

1080 

276 

626 

1440 

376 

678 

1530 

407 

1046 

2160 

628 

1105 

2250 

663 

•••APPROXIMATES  THE  ACTUAL  NUMBER  OF  TRUCKS,  500. 


CUT  NUMBER  V  REDUCED  CAPACITY  ON  ARC  B  TO 

ARC  WILL  BE  RESTOREO  ON  DAY  1. 


D  FROM  1350.  TO 


TRUCKS 

TONS/DAY 

TRUCKS 

IN  INVENTORY 

THROUGHPUT 

ACTIVE 

155 

270 

93 

206 

360 

124 

310 

540 

186 

361 

630 

217 

471 

810 

283 

471**» 

810 

283 

530 

900 

318 

•••APPROXIMATES  THE 

ACTUAL  NUMBER 

OF 

TRUCKS, 

500. 

IT  NUMBER  2  REDUCED 

CAPACITY  ON  ARC 

8  TO 

E 

ARC  WILL  BE  RESTORED 

ON 

DAY  3. 

TRUCKS 

TONS/DAY 

TRUCKS 

IN  INVENTORY 

THROUGHPUT 

ACTIVE 

103 

180 

62 

155 

270 

93 

206 

360 

124 

258 

450 

155 

310 

540 

186 

413 

720 

248 

E  FROM  2700.  TO 


"IU- 


0  *  Y  2 

FLOW  PROFILE  PRIOR  TO  FIRST  CUT 


TRUCKS 

TONS/OAY 

TRUCKS 

IN  INVENTORY 

THROUGHPUT 

ACTIVE 

266 

6)0 

161 

)06 

720 

166 

696*»* 

1170 

299 

336 

1260 

322 

575 

1350 

365 

665 

1690 

5)1 

966 

2070 

593 

•••APPROXIMATES  THE  ACTUAL  NUMB'"  OF  TRUCKS*  900. 


CUT  NUMBER  I  REDUCED  CAPACITY  ON  ARC  B  TO  D  FROM  1330.  TO  0. 

ARC  MILL  BE  RESTORED  ON  DAY  2. 


TRUCKS  TONS/OAY  TRUCKS 

IN  INVENTORY  THROUGHPUT  ACTIVE 

10)  160  62 

133  270  93 

206  360  126 

236  630  133 

310  360  166 

61)  720  26R 


CUT  NUMBER  2  REDUCED  CAPACITY  ON  ARC  A  TO  C  FROM  720.  TO  0. 

ARC  MILL  BE  RESTORED  ON  DAY  2. 

FLOW  TOTALLY  STOPPED  ON  DAY  2 


0  A  Y  3 


FLOM  PROFILE  PRIOR  TO  FIRST  CUT 


TRUCKS 

TONS/OAY 

TRUCKS 

IN  INVENTORY 

THROUGHPUT 

ACTIVE 

306 

720 

166 

365 

BIO 

207 

660«»» 

1080 

2  76 

626 

1660 

376 

676 

1530 

607 

1066 

2160 

628 

1103 

2250 

663 

••APPROXIMATES  THE 

ACTUAL  NUMBER  OF  TRUCKS, 

500. 

CUT  NUMBER  1  REDUCED  CAPACITY  ON  ARC  B  TO  D  FROM  1)50.  TO  0. 

ARC  MILL  BE  RESTORED  ON  DAY  3. 


TRUCKS 

TONS/OAY 

TRUCKS 

IN  INVENTORY 

THROUGHPUT 

ACTIVE 

155 

270 

93 

206 

360 

126 

310 

560 

166 

361 

630 

217 

671 

610 

263 

671«*« 

810 

263 

530 

900 

318 

•••APPROXIMATES  THE  ACTUAL  NUMBER  OF  TRUCKS,  500. 


CUT  NUMBER  2  REDUCED  CAPACITY  ON  ARC  B  TO  E  FROM  2700.  TO  0. 

ARC  MILL  BE  RESTOREO  ON  DAY  3. 


TRUCKS 

TONS/OAY 

TRUCKS 

IN  INVENTORY 

THROUGHPUT 

ACTIVE 

103 

180 

62 

153 

270 

93 

206 

360 

126 

236 

650 

155 

310 

560 

186 

613 

720 

266 

n  n  n  n 


Appendix  B 


INTERDICTION  PROGRAM 


SIBFTC  CIRSEA 

COMMON  /RTRK/  IRFLO 
COMMON  /MARKR  /  mioon 
COMMON  /FACTR/  FACTOR 
COMMON  / R  T I  ME /  INTIME(lOOl) 

COMMON  /NNODE /  I  SINK 
COMMON  /DENOM/  IDENOM 

COMMON  /OR!G/M,N,KB( 12) .NIT .NOT .NN.NP.ll .  I J  .KC .KU .LB . KX , NL  ,  K V 
COMMON  /LFIAG/  K INTRP.KCUT 

COMMON  /FACT/  F1NCOM.  FS TOP . TNSTRK .HRPRDA . TUPRHR 
COMMON  /KDAY/  KPDS »L INK  »L I NK 1 
COMMON  /LBUG/  BUGWRD 

DIMENSION  I  I  (  1001).  I J f 1001  I *KC I  1001  I .KU< 1 001  I • L B (  1 001  I .KX (  1 00 1 ) 
DIMENSION  NLI  501). NP (  5011.NNI  5 01  I ,K V  M.  00 1 ) 

DIMENSION  KCSI  1001  I .KUS<  1  001  )  .L BS ( 1001 ) 

DIMENSION  L0CI2) 

LOGICAL  KINTRP.KCUT 
LOGICAL  BUGWRD 
C 

READ  (5.2)  NTRKS.NPDS.NL  INKS. NPPT  S» IOURCE . I  SI Nk .BuGwRD 
READ  (5.20)  TNSTRK. FINCOM.FSTOP.hRPRDA.TUPRhP 
20  FORMAT  (  3F5.0/2F5.0) 

WRITE  (6.31  NTRKS.NPDS.NLINKS.NPPTS  , IOURCE . I S INK . 

X  TNSTRK.FINCOM.FSTOP.HRPRDA.TUPRHR 

c 

C  INITIALIZATION 

CALL  OKINPT 

C  SAVE  THE  ORIGINAL  DATA 

DO  5000  I  ■  l.N 

KUS ( I  )  -  KU ( I  I 

5000  LBS( I )  ■  LB ( I ) 

C  HERE  FIND  THE  DESTINATION  LINK. 

5001  DO  5002  I  -  l.N 
N2  ■  I  J(  I  ) 

IF  (NN(N2).EQ.  IOuRCE)  GO  TO  50002 

5002  CONTINUE 
50002  ISAVE  •  I 

C  FIND  THE  SUM  OF  THE  DEMAND  RATIOS 

IDENOM  *  0 
DO  5003  I  -  l.N 
N2  «  1 J  I  I  ) 

IF  (  NN(N2).NF.ISINK)  GO  TO  5003 
IDENOM  ■  IDENOM  ♦  LB ( I ) 

5003  CONTINUE 
IDENOM  r  MAXOI 1 , IDENOM) 

MAIN  ROUTINE 

ONCE  FOR  EACH  DAY 

DO  999  KPDS  *  1 . NPDS 
KCUT  «  .FALSE. 

C  THE  NFTWORK  ISNT  CUT  YET 
WRITE  (6.1990)  KPDS 

1990  FORMAT  (  1 H 1  ,  23X.7HD  AY  .12  /  ) 

WRITE  (6.2001) 

2001  FORMAT  (1H  <  1 2X .  32H  FLOW  PROFILE  PRIOR  TO  FIRST  Cul 

C  RESTORE  ALL  REPAIRED  LINKS 
DO  5004  I  ■  l.N 

IF  (  INTI  ME  <  I  1  .LE  .KPDS I  KU(I)  •  KUS(D 

5004  LB ( 1  )  •  0 

NLINK1  •  NLINKS  *  1 
C 

C  NOW  BEGIN  THE  STRIKES  FOR  THIS  CURRENT  DAY 
DO  9990  KLINKS  ■  1.NLINK1 
C  FIND  MAX  FLOW 

KUI ISAVE)  •  KUSI ISAVE) 

LB (I  SAVE)  ■  0 
KC( ISAVE)  *-9999999 

CALL  MNCFIM.N, I  I  .  I J.KC .KU . LB .KX ,NP .NL  .  INFE A S ) 

KC ( I  SAVE )  -  0 

IF  (BUGWRD)  CALL  C JTPijt  (  I  I  .  J  .KC .KU.LB .KX  .NL .NP .  I NFEAS.KV  ) 

C  RESTORE  ALL  TRUE  DEMANDS 
DO  50004  I  «  l.N 
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50004  LSI  I  I  *  *. B S  <  1  ) 

IF  ( ( INFEAS.NE.l > .AND.  « KX ( ISAVE1.NE.0) »  GO  TO  5005 
IF  (  KX I ISAVF ) .FQ.OI  GO  TO  10117 
WRITE  (6.10114) 

10114  FORMAT  (  1H0.51HMIN  DEMANDS  AT  SOME  DESTINATION  CANNOT  BE  SATISFI 
XED  ) 

GO  TO  999 

10117  WRITE  (6.41  KPDS 

004  FORMAT  (  1H  «  6  X . 2  7HF  LOW  TOTALLY  STOPPED  ON  DAY  .  13) 

GO  TO  999 

5005  INC  ■  MAXO( 1 .KX  (  ISAVE  I /NPPTS) 

KSAVE  ■  K X ( I  SAVE ) 

WRITE  (6.2000) 

C  ^E  HAVENT  INTERPOLATED  ON  THf  true  number  of  trucks  yet  for  THE  C0M1N 
C  profile 

KINTRP  •  .TRUE. 

c 

C  FOR  THIS  NETWORK  CONFIGURATION  PLOT  THE  PROFILE 

DO  99990  KPPTS  *  ].NPPTS 
I  FLO  -  I NC*KPP  T S 

C  IEND  IS  TO  TEST  THE  LAST  ITERATION  WITH  ONLY  ONE  WORD 
I  END  ■  NPPTS-KPPTS 

CALL  PROF  I L (  I  SAVE, I  SINK, IFLO.NTRKS. LBS* IDENOM.  IEND, KSAVE I 
IF  ( KCUT )  GO  TO  999 

99990  CONTINUE 

IF  (  .NOT. KINTRP)  WRITE  (6,10134)  NTRKS 
10134  FORMAT  (  48H0  ***APPROX I  MAT ES  THE  ACTUAL  NUMBER  OF  TRUCKS.  ,16. 

X  1H.  ) 

99991  IF  (KLINKS.GT.NLINKS)  GO  TO  9990 
C  NOW  CHOOSE  A  LINK  TO  CUT 

CALL  WOLLMR  (KX.Kij.NN.I  I  ,  I  J,  N  ,KV  .K  L  I NKS  ,LB  .KC  ,M  ,NP  ,NL  .  I  SAVE  ,  LBS  ) 
9990  CONTINUE 
999  CONTINUE 
CALL  EXIT 

C  FORMAT  STATEMENTS 

2  FORMAT  (41 1G,2X.A6,2X,A6.23X,L1 ) 

003  FORMAT ( 10H1  TRUCKS  =  .I6.1X.10H  PERIODS  *  . I  6 . 1  X . 0HL I NK S  «  » 

X  I6.1X.8H  NPPTS  *  .  I6.1X.9HS0URCE  «  , A6 . 1 X , 7HS I NK  «  »A6/ 

X  10H  TONS/ TRK  .F8.4.2X.6HI NCOM  ,  FB. 4 , 2X . 5HST0P  .F8.4.2X, 

X  8HHRS/DAY  .F8.4.2X.6HTU/HR  ,F8.4  ) 

2000  FORMAT  (  1  HO , 6X ,6HT RUCKS . 1 1 X .8HT0NS /DA Y , 1 7X .6MTRUCKS/ 

X  4X.12HIN  INVENTORY. 7X  ,1 OHTHROUGHPUT ,  1 6X , 6HACT I VE  ) 

END 

SIBFTC  WOLLMR 

SUBROUTINE  WOLLMR ( KY ,KU .NN » I  I , I J ,N ,KV ,KL I NKS »LB ,  KC  ,M ,NP , NL , I  SAVE , 
X  LBS) 

COMMON  /RTRK/  IRELO 
COMMON  /NNODE /  I  SINK 

COMMON  /DENOM/  IDENOM 

COMMON  /FACTR/  FACTOR 

COMMON  /FACT/  F  I  NCOM,  FSTOP » TNSTRK ,HRPRDA , TUPRHR 

COMMON  /MARKR /  I T  < 1001 ) 

COMMON  /LBUG/  BuGWRD 
COMMON  /RT IMF/  I  NT IME 1 I0OJ  > 

COMMON  /KDAY /  KPDS ,L I NK . L  I  NK  1 
DIMENSION  LB ( 1  I  ,KC  ( 1  )  .LBS ( 1 ) 

DIMENSION  KY( 1 ) «KU ( 1  )  «NN(l),t I ( 1  I • 1 J  (  1  ) , K V ( 1 ) 

DIMENSION  LOC ( 2  I 
LOGICAL  BUGWRD 
WRITE  (6,10017) 

10017  FORMAT  (  1H0) 

LINK  ■  0 

MFLO  ■  999999 
ITRO  •  0 
LB ( I  SAVE  )  ■  0 
KC ( l SAVE )  ■  -9999999 
KU ( I  SAVE )  ■  IRFLO 

I F ( IRFLO.EO.O)  KU ( I  SAVE  I  ■  9999999 
DO  1  I  -  1  »N 

I  T ( I  1  ■  M I  NO ( K V ( I  I , KU ( I ) ) 

N2  ■  I J( I ) 

IF(NN(N2).NE.ISINK>  go  TO  1 
KLB  ■  L8S< 1 )*!RFL0/I DENOM 
LB ( I )  •  Ml  NO ( LBS (  I  1  «KLB ) 

LB ( I  SAVE )  «  LB ( I  SAVE )  ♦  L8III 
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001  CONTINUE 

CALL  MNCF  (M.N.  I  I .  I  J.KC  .KIJ.LB  .ICY  .NP*NL  •  INFEAS) 

OO  2  IMAX  -  l.N 

004  IE  (  III IMAXt.EQ.O)  GO  TO  ? 

I  EX  ■  0 
00  3  I  •  l.N 

C  MARK  OUT  OF  CONSIDERATION  ALL  ARCS  WITH  ZERO  FLOW 
IE  (  KYI  I ) • EQ.  01  IT(I)  *  0 
003  I E X  ■  I  EX  ♦  I  T ( I ) 

C 

IE  <  IEX.EQ.0)  GO  TO  35 
I  TEMP  •  KU I IMAX ) 

KU ( IMAX)  *  0 
KC< ISAVE  )  -  -9999999 

C  FIND  MAX  FLOW  IN  NETWORK  WTTH  LINK  IMAX  OUT 

CALL  MNCF  I  M.N.  I  1  .IJ.KC.KU.L8.KY.NP.NL.INEi:  AS) 

KC ( I  SAVE )  *  0 
KUIIMAX)  *  I  TEMP 
IE  (  INEEAS.NE.1 !  GO  TO  6 
LINK  ■  IMAX 
GO  TO  35 
C 

C  IE  RESULTING  flow  IS  GREATER  than  TEST  FLOW  SKIP  THIS  LINK,  IF  LESS  KEEP  THE 
C  LINK.  AND  IF  EQUAL  F I  NO  OUT  THE  TOTAL  NUMBER  OF  TRUCKS  REQUIRED. 


IF 

THEY 

ARE  LESS  FORGFT  IT. 

006 

IF  ( 

K Y (  I  SAVE) .GT. MFLO) 

GO 

TO 

2 

CALL 

INPRD(KY.KC.N.NUTRQ) 

IF  ( 

KYI  ISAVE). LT. MFLO) 

GO 

TO 

7 

IFINUTRQ.LT.  1TRGI  GO  TO 

2 

007 

LINK 

■  IMAX 

MFLO 

■  K Y (  I  SAVE  1 

I  TRQ 

»  NUTRQ 

002 

CONTINUE 

035 

ZK 

•  KU(LINK) 

IF  (LINK.NE.O)  GO  TO  10116 
WRITE  (6.10117) 

10117  FORMAT ( 23H  NO  LINKS  ARE  VULNERABLE.) 

RETURN 

1C116  KU( LINK )  «  0 
AK  ■  KVILINKI 

KOUT  «  AK/ ( TUPRHR*HRPRDA ) w.5 
INTIME(LINK)  ■  KPDS  ♦  KOUT 
A  «  ZK/FACTOR 
OK  «  KU(LINK) 

B  »  QK/FACTOR 
N1  «  IKLINK) 

N2  •  IJUINKI 

WRITE  (6.9)  KLINKS.NN(Nl  I  ,NN(N2I . A .8 . 1  NT  I  ME ( LINK ) 

009  FORMAT ( 1 2H  CUT  NUMBER  .I2.23H  REDUCED  CAPACITY  ON  ARC  .A6.4H  TO  . 

XA6.  6H  FROM  .F7.0.4H  TO  . F 7.  0 / 1 5 X . 2$H  ARC  WILL  BE  RESTORED  ON  DAY 
X.I3.1H.) 

C  THE  NEXT  SECTION  OF  CODE  IS  TO  TAKE  OUT  LINKS  IN  THE  OPPOSITE  DIRECTION  IF 
C  THEY  EXIST 

IRFLO  -  0 

IF  (  IEX.EQ. 999999)  GO  TO  365 
DO  36  I  -  l.N 
Ml  *  II(I) 

M2  ■  IJU) 

I F (  NNIM1 ) .EQ.NNI N? ) .AND.  NN ( M2 ) . EQ .NN  I N]  I.AND.KVI I ) . EO.K V ( L I NK ) ) 

X  GO  TO  37 
036  CONTINUE 
RETURN 

365  CONTINUE 

LINK  «  ITTMP 
RETURN 

037  ITTMP  »  LINK 
LINK  »  I 
LINK1  «  I 
I  EX  ■  999999 
GO  TO  35 

END 


C 
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*  1 BF TC  PROF IL 

SUBROUTINE  PROFILI  l  SAVE ,  I S  l  NIC  ,  IFlO.NTRKS  ,LBS.  IDENOM,  I  END.KSAvE  I 
COMMON  /LFLAG/  ICtNTRP.KCUT 

COMMON  /OR  I  G/M.N. KB (  12)  .NIT.NOT  .NN  ,NP  «  1  I  •  1  J.ICC  .ICUtLB  *«CX  »NL  .K  V 
COMMON  /LBUG/  BUGWRD 

01  MENS  ION  I  I  <  1001  )  .  I  J(  1001 1  .KC(  1001  )  .KUU001  I  .LB(  1001 )  .KX(  1001  ) 
DIMENSION  NL(  S01).NP(  501). NN(  501 ) »K V ( 1 00 1 ) 

DIMENSION  LBS ( 1001 ) 

DIMENSION  LOC I  2 ) 

LOGICAL  KtNTRP.KCUT 
LOGICAL  BUGWRD 
JT  ■  1 

010  LBIISAVE)  ■  0 

DO  1000  I  ■  l.N 
N2  >  I  Jill 

IF  I  NNt N2 ) .NE. I S  INK )  GO  TO  1000 
KLB  «  LBSI I I*1FL0/1 OENOM 
LB( 1 )  -  MINOILBSI  I  )  • KLB ) 

1000  CONTINUE 

LBIISAVE)  -  IFLO 

CALL  MNCF (M.N.II.IJ.KC.KU. LB. KX.NP.NL. INFERS) 

IF  I  INFERS. NE.l  )  GO  TO  1010 
WRITE  (6.2010) 

2010  FORMAT  (  1  HO .51 HM I N  DEMANDS  CANNOT  BE  SATISFIED  AT  SOME  DESTINATI 

XON  1 
IF  (BUGWRD) 

XCALL  OUTPUT (II.  I  J ,KC . Ku .LB  .  KX , NL .NP . I NFE AS .K V ) 

KCUT  -.TRUE. 

RETURN 

1010  CALL  INPRDI KX.KC  .N.  I  TRQ) 

IF  (BUGWRD)  Call  OUTPUT! II . IJ. KC.KU.LB.KX.NL. NP. INFERS. KV) 

CALL  JADJ ( I TRO. IFL0.JTRQ.JFL0) 

IF  ( ( JTRQ.GE.NTRKS).AND.(KINTRP) ) 

X  CALL  INTRPI JTRO.NTRKS.IFLO. 

X  LTRQ.LFLO. I  SAVE, IDENOM.LBS.ISINK) 

L  TRQ  «  JTRQ 
L FLO  •  I  FLO 

WRITE  (6.2000)  JTRQ.JFLO.ITRQ 
2000  FORMAT  (  1 H  .  5X ,  I6.10X.18.15X.I9) 

IF  (  IEND.NE •  0)  GO  TO  1020 
IF  ( I  FLO. EQ.KSAVE I  GO  TO  1015 
IFLO  ■  K SAVE 
GO  TO  10 
1015  LFLO  -  0 
LTRQ  »  0 
RETURN 

1020  IF  (  JT.EQ.  1)  GO  TO  1030 
RETURN 
1030  JT  »  2 

IFLO  •  IFLO  ♦  1 

GO  TO  10 

END 

tIBFTC  INTRP 

SUBROUTINE  INTRPI  I  TRO,  NTRK S . ! FLO . L TRO .LFLO  .  I  SA VE . I DENOM ,LBS .  I  S 1 NK ) 
COMMON  /RTRK/  IRFLO 
COMMON  /LBUG/  BUGWRD 
COMMON  /LFLAG/  KINTRP.KCUT 

COMMON  /ORIG/M,N,KB( 12 ) .N I T .NOT .NN.NP. I  I  .  I  J.KC .KU.LB »KX .NL »KV 
DIMENSION  I  I  (  1001  )  .  I  J(  1001  )  ,KC(  10011  .KUU001 )  .LBl  1001  )  .KX(  1001  ) 
DIMENSION  NL (  SOI ) ,NP<  501). NNt  501 ) »KV( 1001 ) 

DIMENSION  KCS(  1001)  ,KUS( 1001) ,lBS( 1001 1 
DIMENSION  LOC  <  2  ) 

LOGICAL  KINTRP.KCUT 

LOGICAL  BUGWRD 

INTEGER  OLOFLO 

IF  (  I TRQ.EQ.NTRKS  I  GO  TO  2003 

IRFLO  «  0 

NUFLO  «  LFLO*(1FLO-LFLO)»(NTRKS-LTRQ>/I!TRO-LTAQ) 

IF  (  NUFLO. EO. LFLO)  NUFLO  «  NUFLO  ♦  1 
500  DO  1000  I  ■  l.N 
N2  *  I J  <  I  I 

IF  (  NN(N?).NE.IStNK)  GO  TO  1000 
KLB  >  LBSI I ) *NUFLO/ IOENOM 
LB ( I  I  «  MINO(KLB»LBS( I )  1 
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1000  CONTINUE 

NUT  NIC  ■  0 

LB ( I  SAVE )  •  NUFLO 

CALL  MNCF (M.N. I | .IJ.KC.KU.LB.KX.NP.NL. 1NFEAS I 

CALL  IMPRDUX.KC.N.NUTRKI 

CALL  JAOJINUTRK  .NUFLO.JTRQ .JFLO) 

IRFLO  •  NUFLO 

IF  (  JTRQ.GE.NTRKS-  501  GO  TO  2003 
WRITE  (6*2001  I  JTRO.JFLO.NUTRK 
2001  FORMAT  (  1H  ,5X , 1 6 • 10X , 1 8 . 15X  .  1 9  I 
NUTEM1  •  NUFLO  ♦  ( I FLO-NUFLO ) /3 
NUTEM2  ■  NUFLO  ♦  1 
NUFLO  ■  NAXO (NUTEM1 *NUT£M2  ) 

GO  TO  300 

2003  WRITE  (  6.20001  JTRQ* JFLO* NUT RK 

2000  FORMAT  (  1H  ,  5X,  I6*3Ha**.7X*I8.15X»l9) 

KlNTRP-  .FALSE. 

RETURN 

END 

SIBFTC  OKINPT 

SUBROUTINE  OKINPT 
COMMON  /FACTR/  FACTOR 

COMMON  /FACT/  F I NCOM. F STOP  .  TNSTRK .HRPRDA.TUPRHR 
COMMON  /ORIG/M.NtKBl  12)  .N  IT  .NOT  .NN.NP.  1 1  .  I  J.KC  .KU.LB.KX  .NL  ,KV 
DIMENSION  NL<  301)*NP(  301I.NNI  501) 

DIMENSION  1 1 ( 1001 ) • I J (  1001  I .KC< 1001  )*KU( 1001 >  *LB(  1001 ) *KX( 1001 1 
DIMENSION  KA(12).KV(1001) 

C  INPUT  TAPE 
NIT  «  3 

C  OUTPUT  TAPE 
NOT  »  6 

C  MAXIMUM  ARCS 

MAXN  •  1000 
C  MAXIMUM  NODES 
MAXM  a  S00 

FACTOR  a  1.0/(TNSTRX*TUPRHR*HRPRDAI 
C  PREPARE  TO  READ  DATA 

100  READ  (NIT. 90  I  ( KA ( I ) a  I *1 1 1 2  ) 

CALL  NUM( K A* 36HPAUSE  SAVE  READY  ARCS  PUNCHAPUNCMN.KKK  ) 

IF  ( KKK . EQ»  2  .OR.  KKK.EQ.3I  GO  TO  102 
WR I TE ( NOT .91 J  (KA ( I) • I al , 1 2 | 

IF  t  KKK.GT .6 )  GO  TO  100 
IF  (KKK.EQ.A)  GO  TO  200 
IF  (KICK.EQ.5l  GO  TO  160 
IF  (KKK.EQ.6)  GO  TO  170 
160  CALL  EXIT 
C  PUNCHING 

160  PUNCH  79. IKAI I | .|al  .11  ) 

DO  161  J  *  l.N 

N 1  ■  I  I ( J) 

N2  •  IJ(JI 

PUNCH  78.NNIN1 ) .NNIN2) .KC( J) *KU (  J )  .LB ( J ) *KX ( J ) 

161  CONTINUE 
PUNCH  T7 
GO  TO  100 

170  PUNCH  76 

DO  171  I  ■  l.M 

PUNCH  7S.NNI I  I ,NP< I  I 

171  CONTINUE 
PUNCH  7T 
GO  TO  100 

102  WRITE  (NOT .97 )  (KA( 1 1 . lai,  12) 

IF  (RKK.E0.2)  GO  TO  101 
C  ZERO  NOOE  PRICES 

110  DO  111  I ■  1  .MAXM 

NPt I )  a  0 

111  CONTINUE 

C  READ  TITLE  CARD 

101  READ  (NIT. 90)  (KB( I ) . I»1 .12) 

WRITE  (NOT. 91)  ( KBI I ) » I  a] , 12 > 

IF  (KKK.E0.2)  GO  TO  501 

GO  TO  100 
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C  HEAD  ARCS 

200  M  «  0 

WRITE  ( 6  •  1 0 1 26 ) 

10126  FORMAT (5H0ARCS16X.4HC0ST6X.  5HUPPE  R6X «  5HL0WER •  10X.4HTIME/ 

X  5H  FROM .3X.2HT0.46X.9HTO  REPAIR/1XI 
00  201  1-1,10000 
*3?  FORMAT  (‘JA6.2X.4llO.131 

REACINI  T  .<32  I  KD.III  IltIJI  ll.KCI  I  )  ».TU (  I  )  •  JLBT  .XX t I  I .XVt I ) 

CALL  NUMIXD.  24HFND  NODES  COMPUT  .XXX) 

IF  (XXX. NE. 4  )  GO  TO  209 
ALBT  •  JLBT 
LB ( I )  «  ALBT»FACTOR 

WRITE  (6. 10136) III  I) . I J< I )  .XC 111  .XUI  II  .LBI I ) ,  XVII) 

10135  F ORMAT (2A7.3I11.  11X.I3) 

201  CONTINUE 

209  WRITE  I  NOT  .931  XD. I  I  I  I ) . I J I  1 1 

?K  N  «  1-1 

220  IFIXXK.NE.5)  GO  TO  TOO 

299  WRITE  I  NOT  »  96 ) 

CALL  EXIT 

C  NUMBER  THE  NODES 

300  DO  301  L  ■ 1 *N 

NN(M-M)  «  I  1  I L ) 

I  I  I  L I  -  NODENO I  I  I  I L )  ) 

M  *  MAXOtM.lt (LI  I 

301  CONTINUE 

DO  302  L-l.N 

NNIM41)  «  IJIL) 

I J I L I  ■  NODFNOI IJIL) I 
M  •  MAXOIM, IJIL '  I 

302  CONTINUE 

C  SXIP  NODE  READING  IF  NODE  CONTROL  WASN»T  READ 
IF  (XXX. NE. 2)  GO  TO  500 
C  NODF  READING 

400  READ  (NIT. 81)  XD.KE.NPX 

CALL  NUMIXD.  18HEND  COMPUT, XXX) 

IF  (XXX.NE.2I  GO  TO  410 
N 1  «  NOOENOIXE) 

IF  I  Nl.GT.M  )  GO  TO  400 

420  NPINl)  ■  NPX 
GO  TO  400 

410  WRITE  (NOT. 82)  XP.XE 

IF  (XXX. EQ. 4)  GO  TO  299 
GO  TO  500 
C  REAOER 

500  IF  (XXX. EQ. 3)  GO  TO  600 

501  READ  (NIT. 83)  < XA < l  ) . I -1 , 8  I 

CALL  NUMI XA.24HALTER  COMPUT  NODES  .XXX) 

IF  ( XXX. EQ. 1  .OR.  XXX. EQ. 3  I  GO  TO  515 
WRITE  (NOT. 93)  XA ( 1  I .X A ( 2 )  . XA <  3 ) 

560  IF  (XXX.EQ.2l  GO  TO  600 
IF  (XXX. EQ. 5)  GO  TO  501 
IF  (XXX. EQ. 4)  GO  TO  400 

515  WRITE  (NOT. 841  ( X  A  C  1  I , I ■ 1  ,  8  ) 

X  A ( 4 )  •  MAX0IXAI4) .1  I 
N1  •  NOOENO ( X A ( 2 ) ) 

NNINl  )  -  X A ( 2  I 
M  ■  MAXOIM.Nl ) 

N2  «  NOOENO (XA(  31) 

NNIN2)  ■  X A  (  3 ) 

M  ■  MAXOIK.N2) 

DO  520  I  -l.N 

IF  (Nl.NE.II(I).OR.N2*NE.IJ< I)  )  GO  TO  520 

X  A  (  4  |  •  X  A  (  4  )  -  1 

IF  ( XA ( 4  )  . EQ, 0 )  GO  TO  530 

520  CONTINUE 

WRITE  (NOT, 85) 

N  «  N*1 
I  -  N 

530  HID  ■  N1 
I J ( I )  ■  N2 
KC ( I )  -  X A ( 5  I 
XU ( I >  ■  XA (6) 

LBI I )  ■  X A I  7 1 

XXI  I )  *  X A ( 8 1  ♦  XX( I | 


rv  r\ 
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go  TO  501 

600  WRITE  (NOT. 94 1  N  »M 
C  TESTS 

IF  (N.LE.MAXN)  GO  TO  601 
WRITE  (NOT .95 ) 

CALL  EXIT 

601  IF  (N.LE.HAXH)  GO  TO  320 
WRITE  (NOT .98  I 

CALL  EXIT 

C  LOOK  FOR  DEAD  NODES 
320  00  303  !•!  «M 


DO  304  L  ■ 

1  .N 

IF(IKL) 

.EQ. 

>  I  ) 

GO 

TO 

305 

304 

continue 

WRITE (NOT. 

99) 

NNI 

1  I  ) 

305 

DO  306  L  ■ 

1  .N 

IFIIJ(L) 

.EO. 

.11 

GO 

TO 

303 

306 

continue 

WRITE (NOT. 

80) 

NNI 

:  I  ) 

303  CONTINUE 

DO  611  I •  1  «M 
NL ( I >  •  0 

611  CONTINUE 

C  CALCULATE  flows 
DO  612  J»1.N 

IF  (RU(JI»GE.L8( J)  I  GO  TO  609 
LUP  •  KU  (  J ) 

KU(J)  ■  LB ( J I 
LB(J)  ■  LUP 
609  N1  ■  I  I  (  J ) 

N2  «  IJIJI 

NL(Nl)  •  NL(N1)  -  KXIJI 
NLIN2)  ■  NL ( N2 )  ♦  KXIJI 

612  CONTINUE 

DO  603  l-l.M 

IF  (NL(I).NE.Ol  WRITE  (NOT. 661  NN(I),NL(I( 

605  CONTINUE 
C  DO  THE  ALGORITHM 
RETURN 

75  FORMAT (6X.A6. 1181 

76  FORMAT (6HN00ES  I 

77  FORMAT (6HEND  I 

78  FORMATI6X.2A6.2X. 4110) 

79  FORMAT (6HARCS  .11A6I 

80  FORMAT! 24H  NO  ARC  ENDS  AT  NODE  A6 I 

81  FORMAT ( 2A6 .118) 

82  FORMAT! IX. 2A6 ) 

83  FORMATI3A6. 12.41 10) 

84  FORMAT! IX.3A6. 14.41101 

85  FORMAT ( 26H  AUGMENT  ARCS  BY  ABOVE  ARC) 

86  FORMAT ( 6H  NODE  .A6.28H  NON-CONSERVATIVE.  NET  FLOW-114) 

90  FORMAT ( 12A6  I 

91  FORMAT ( 1H1 . 12A6  ) 

93  FORMAT (IX. 3A6 1 

94  FORMAT ( 18H  THIS  PROBLEM  HAS  I12.9H  ARCS  AND. 112. 6H  NODES  ) 

95  FORMAT ( 14H  TOO  MANY  ARCS) 

96  FORMAT ( 34H  ILLEGAL  OATA  CARD  OR  CONTROL  CARD) 

97  FORMAT(1H0.12A6) 

98  FORMAT (15H  TOO  MANY  NODES) 

99  FORMAT  ( 24H  NO  ARC  BEGINS  AT  NODE  A6 ) 

END 

HBFTC  MNCF 

SUBROUTINE  MNCF (NODES. ARCS.  I  .J.  COST. MI jLO.FLOW.PI . NA. JNFEAS) 
INTEGER  NODES. ARCS. I. J. COST. HI  .LO. FLOW. PI .NA.INFEAS 
DIMENSION  I ( 3000) .J( 3000) .COST (  3000  I .Ml  (  3000 ) »L0(  3000 > .FLOW! 3000 ) 
DIMENSION  LOC ( 2 ) 

DIMENSION  PM  1000) #NA< 10001 
C 
C 

DEFINITION  OF  CALLING  SEQUENCE 


* 


C 

c 


NAME  USE 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


NODES  NUMBER  OF  NOOES 

ARCS  NUMBER  OF  ARCS 

I  LIST  OF  FIRST  NOOES 

J  UST  OF  SECOND  NOOES 

COST  UNIT  COST  OF  FLOW  ON  ARCS 
HI  UPRER  BOUNDS  FOR  ARCS 

LO  LOWER  BOUNDS  FOR  ARCS 

flow  amount  of  flow  in  arcs 
PI  NOOE  PRICES 

NA  SCRATCH  list  FOR  NOOE  LABELING 

INFEAS  Flag  denoting  the  condition  of  output 

BEGIN 

INTEGER  A.AA.N.SRC.SNK.DEL  < INF  ,C  .  AOK ,COK *N 1 ,N2 • INC .L ABEL 
INF  •  34JS9738367 
AOK  ■  0 


C  LOOK  FOR  AN  OUT  OF  KILTER  ARC 
DO  90  AA*  1 .  ARCS 
100  N1  *  I ( AA ) 

N2  ■  J I AA  I 

C  »  COST(AA)  ♦PI (Nil  -  P1(N2) 

C  EXIT  IF  SENSE  SWITCH  5  IS  DOWN 
30  CALL  SSWTCH(3. LABEL) 

IF  (LABEL. EO. 2)  GO  TO  40 

INFEAS  ■  2 

RETURN 

40  IF(FLOW(AA).LT.LO(AA|.OR. (C.LT.O.ANO.FLOWIAA) .LT.Hl I AA) )  ) GOTO  50 
IF(FLOWC AAl.GT.MI (AA).OR. ( C . GT  .  0. AND . FLOW  (  A A 1 .GT . LO ( AA ) )  IGOTO  60 
90  CONTINUE 

C  NO  OUT  OF  KILTER  ARCS  LEFT 

INFEAS  ■  0 
RETURN 

C  OUT  OF  KILTER  ARC  FOUND 
50  SRC  ■  J ( AA ) 

SNK  «  I ( A A  ( 

LABEL  ■  >A A 
GO  TO  200 
60  SRC  ■  I ( AA I 
SNK  ■  J ( AA ) 

LABEL  •  -AA 

C  SAVE  LABELS  IF  LAST  OPERATION  WAS  INCREASING  NOOE  PRICES  ON  THIS  ARC 

200  IF  (  AA.EO.AOK.AND.NAISRC ) .NE.O)  GO  TO  203 
DO  201  N  •  1. NODES 

NA ( N )  «  0 

201  CONTINUE 
AOK  »AA 

205  COK  «  C 

NA(SRC)  ■  LABEL 
C  LABEL 

210  LABEL  •  0 

DO  250  A  « 1 .ARCS 


N 1  •  1(A) 

IF  (Nl.LT.Ol  GO  TO  250 
N2  »  J ( A ) 

IF  (NA(N1).EO.O.ANO.NA(N2).EO.O)  GO  TO  250 
IF  (NA(Nl). NE.O. AND. NAIN2 ). NE.O)  GO  TO  245 
C  ■»  COST ( A )  ♦  PI (Nil  -  PI ( N 2  ) 

IF  <  NA ( N 1 l.EO.Ol  GO  TO  220 

I F ( FLOW ( A>.GE.HI(A).OR.( FLOW ( A  I «GE.LO( A ) . AND.C.GT.O ) )  GO  TO  245 
NA ( N2 I  -  A 
GO  TO  240 

220  IF(FLOW(A).LE.LO(A).OR.(FLOW(A).LE.HI(A).ANO*C.LT.O) )  GO  TO  243 
NA(Nl)  •  -A 
240  LABEL  »  1 

C  NODE  LA8ELE0.  TEST  FOR  BREAKTHRU 
IF  (NA(SNK).NE.O)  GO  TO  260 
245  I ( A)  ■  -N1 
250  CONTINUE 

C  GO  BACK  AND  DO  MORE  LABELING  IF  SOME  NOOE  WAS  LABELED  ON  LAST  PASS 
IF  (LABEL. NE.O)  GO  TO  210 
C  RESTORE  POSITIVE  SIGNS  TO  FIRST  NOOE  LIST 
260  DO  270  A  »  l.ARCS 

I (A)  >  I ABS (1(A)  I 
270  CONTINUE 


r»  n  r\ 
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C  IF  NO  LABELING  DONE  UN  LAST  PASS.  GO  TO  INCEASE  PIE 
IF  (LABEL. EU.O)  GO  TO  400 
C  BR  r AK  THRU  »  FIND  THF  INCREMENT 
300  INC  ■  INF 

C  FOLLOW  PATH  BACK  FROM  SOURCE 
N  «  SRC 
3 1 0  A  *  I A8S ( NA ( N  )  | 

IF  (NAIN) .LT.O)  GO  TO  315 
N 2  «  1(A) 

C  *  COST!  A )  -  PKN)  ♦  PKN2I 

IF  (C.GT.O)  INC  *  M I  NO  < INC.LO(A)-FLOWIA)  ) 

IF  (C.LE.OI  INC  =  MINO( INC.HI  ( AI-FLOW(A)  ) 

GO  TO  340 
3 1 S  N2  *  J ( A  ) 

C  =  COST ( A  I 
3 ?0  IF  (C.LT.OI 
IF  (C.GF.O) 

340  N  =  N? 

IF  (N.NF.SRC)  GO  TO  310 
C  INCRFMENT  ARCS 

3S0  A  =  I ABS ( NA ( N  )  ) 

IF  (NA(N).LT.O)  GO  TO  360 
F  LOW ( A )  «  FLOW ( A  )  +  INC 


+  PI (N)  -  PI ( N2 ) 

INC  =  MINOI INC«FLOw<A|-HJ  (A)  ) 

INC  =  MINO ( INC .FLOWI A )-LO( A )  ) 


N  *  1(A) 

GO  TO  370 

360  FLOW ( A )  «  FLOW ( A  )  -  INC 
N  ■  J  (  A  ) 

370  IF  (  N.NE.SRC  )  GO  TO  350 
C  FLOW  INCREMENTED.  RFTljRN  TO  KILTER  TEST 
NAIN)  «  0 
GO  TO  100 
C  CHANGE  PI 
400  DEL  «  INF 
C  FIND  INCREMENT 

DO  420  A* 1 .ARCS 
N 1  «  1(A) 

N2  ■  J  (  A  ) 

IF  (NA(Nl).EQ.O  .AND.  NA(N2).EQ.O)  GO  TO  420 
IF  (NA(Nl).NE.O  .AND.  NA(N2).NE.O)  GO  TO  420 
C  *  COST ( A )  ♦  PI (N1 )  -  PIIN2I 

IF  <NA(N2I.EQ.O  .AND.  FLOW ( A ) .L T .H l < A )  )  DEL-M I  NO ( DEL .C ) 

IF  (NAIN2I.NF.0  .AND.  FLOW! A  I .GT.LO( A )  )  DEL«MI NO ( DEL  «-C ) 

420  CONTINUE 

IF  (DEL. NE. INF )  GO  TO  430 

IF  (FLOWI  AA)  .EO.LO(AA)  .OR.  FLOW  (  AA  )  .EG)  .H  I  (  A  A  )  )  GO  TO  425 
C  INFEASIBLE  SOLUTION 
INFEAS  «  1 
RETURN 

C  INCREASE  PI 

425  DEL  *  IABS(COK) 

430  DO  450  N  »1 .NODES 

IF  (NA(N).EO.O)  PI(N)  «  PI(N)  ♦  DEL 
450  CONTINUE 

C  GO  BACK  TO  KILTER  TEST 
GO  TO  100 
END 

5IBFTC  OUTPUT 

SUBROUTINE  OUTPUT  (  I  I  *  I  J*KC  . Kl) *LB * KX «NL «NP »  I NFE AS »K V ) 

COMMON  /RTIME/  I  NT  I  ME ( 1001  ) 

COMMON  /LBUG/  BUGWRD 
COMMON  /MARKR  t  IT (1001) 

COMMON  /ORIG/M.N.KBI 12) .NIT .NOT .NN 
DIMENSION  NL(  1001  )  .NPU001  )  .NN(  1001  ) 

DIMENSION  I  I (  3001  )  .  I J  (  3001  ) .KCI3001 ) .KU< 300 1  1  »LB ( 3001  I »KX( 3001 ) 
DIMENSION  KV( 1001  ) 

LOGICAL  BUGWRD 

CARCS  COST  UPPER  LOWER  X  CBAR 

CAAAAAA  AAAAAA  ...I  ...I  ...I  ...I  ...I 


KILTER  NUMBER 
I 

92  FORMAT (5H0ARCS16X.4HC0ST6X.5HUPPER5X .5HL0WER  7X.1HX  8X .4HCBAR3X > 

1  13HKILTER  NUMBER. 5X. 4HVULN, 4X. 6HINT IME.2X.2HIT/1X) 

93  F0RMATI2A7.4I 10, 18  .  I  7. 5X  ,  I 3.4X , I 6.2X , I  3  I 
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098  F0RMAT(2A7,4l 10. 18  .  1  7 , 3X  .  I  3 . 4X  .  1 6 . 2X  .  1 3  1 3HCUT  ) 

WRITE  (NOT  *92  I 
NOUT  *  0 
KILTER  .  0 
MA  ■  0 
MB  •  0 

DO  901  J-l.N 
N1  •  II (J) 

N2  •  IJ(J) 

KCK  •  KC(J!  ♦  NP(N1)  -  NPIN2! 

IF  (KCK.GT.O)  KIL  ■  I ABS ( KX ( J I -LB UII 
IF  (KCK.LT.O)  x I L  -  IABS(KX(J)-KU(Jt I 

IF  IKCK.EQ.O)  KIL  -  MAXO < MAXO ( 0 .KX IJ I -KU ( Jit . LB ( J » -KX ( J ) ) 

IF  (KIL.NE.O)  NOUT  ■  NOUT  *  1 

KILTER  ■  KILTER  ♦  KIL 

CALL  MADOf  MA .MB  *KC ( J I .KX ( J) ) 

IF  ( INFEAS.EO.l  .AND. (NL ( N 1  I *NL ( N2 ) > . EO.O. AND. ( NL ( N 1 ) ♦NL ( N2 ) I 
C  .NE.O)  GO  TO  910 
IF  (KX(J).GT.O) 

XWR I TE (NOT .93 )  NN(Nl) »NN(N2)*KC(J) »KU(  J I *LB ( J ) . KX ( J I 

X  .KCK.KIL.KV(J).INT1ME(J).IT(J) 

GO  TO  901 

910  WRITE  I  NOT .98 )  NN ( N 1) .NN ( N2 I »KC ( J I »KU( J I .LB ( J ) »KX ( J ) 

X  »KCK*KIL»KV(J|»INTIME(J)tlT(J) 

901  CONTINUE 
ME  »  1E*09 
MD  ■  MB/ME 

MB  ■  MB  -MD*ME 
Ml  »  MA  ♦  MD 

IF  (  Ml. EO.O. OR. (  ISI GN (1 .Ml  ) . EQ< I  SI GN( 1 *MB ) )  )  GO  TO  922 
MB  •  MB  4  I S I GN ( ME . Ml  I 
Ml  -  Ml  -  ISlGN(l.Ml) 

922  IF  ( Ml .EO.O )  Ml  ■  IS1GNIM1.MBS 
M2  •  MB  >  I  SIGN  (  ME • MB  I 
CALL  INPRD<  KX.KC.N.MT  ) 

WRITE (NOT. 94)  NOUT . K I L TER  .MT .Ml ,M2 
94  FORMATI1HO.I4.41H  ARCS  OUT  OF  KILTER.  TOTAL  KILTER  NUMBER-18. 
X2X . 1 7HTRUCKS  REQUIRED-  .3(16.1X11 
RETURN 

WRITE  (NOT. 93) 

093  FORMAT  (  1H1.12H  NODE  PRICES/1X) 

MINI  -  NP (  1  ) 

DO  940  I  >  2.M 

MINI  ■  MINO ( NP ( I ) .MINI  I 
940  CONTINUE 

DO  902  I-l.M 

NP( I )  «  I ABS ( NP ( I  I  -  MINI  ) 

IF  (  INFEAS.HF.l  .OR.  NHIJ.EQ.O  )  GO  TO  903 
WRITE  (NOT, 991  NN ( I ) .NP ( I  ) 

99  FORMAT! A7. I 13.8H  LABELED) 

GO  TO  902 

903  WRITE  (NOT. 96)  NN(l).NP(!l 
96  FORMAT ( A7. 1 13 ) 

902  CONTINUE 
RETURN 
END 

HBFTC  NOOENO 

FUNCTION  NOOENO  (III) 

COMMON  /0R1G/M.N.KBI 1 2 )  .N  I T  .NOT .NN( 1000 ) 

NOOENO  ■  M  4  1 
IF  (M.CO.O)  RETURN 
DO  10C1  I-l.M 

IF  (III .EO.NN( I )  )  NOOENO  ■  I 
1001  CONTINUE 
RETURN 
ENO 

SIBFTC  NUM 

SUBROUTINE  NUM  (KA.KB.LI 
DIMENSION  KB( 1) 

00  1  1  •  1.100 

I f  (KB( I ) »EQ. -34339738367  )  GO  TO  2 
IF (KA.EQ.KBt I ) )  GO  TO  2 

1  CONTINUE 
I  ■  1 

2  L-I 
RETURN 
END 
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tlBMAP  MAOD 

ENTRY  MAOD 

MADD  SAVE  4 

LAC  MADD. 4 

LDQ*  4.4 

MP  Y •  5,4 

STO  MA 

STQ  MB 

CLA  MB 

ADO*  3,4 

STO*  3,4 

ARS  35 

ADD  MA 

ADD*  2,4 

STO*  2.4 

RETURN  MADD 

MA  BSS  1 

MB  BSS  1 

END 

1IBFTC  INPRD 

SUBROUTINE  I NPRO ( MZ A . NZ A . N ,  I  TEMP  I 
DIMENSION  M2  A ( 1  |  ,NZA( 1 ) 

I  TEMP  •  0 
DO  1  I  •  1  ,N 

001  ITEMP  »  ITEMP  ♦  M2  A ( I | *N  Z  A ( I  | 

RETURN 

END 

MBFTC  JADJ 

SUBROUTINE  JAD J ( I T RO . 1  FLO . J TRQ , JF LO ) 

COMMON  /FACT/  FINCOM,  F S TOP , TN ST RK , HPPRDA • TUPRHR 
A  «  I TRQ 

ATRUCK  ■  A*FST OP /FINCOM 

C  ATRUCK  IS  THF  ACTUAL  NUMBER  OF  TRUCKS  REQUIRED 
B  *  IFLO 

BFLO  *  B*TUPRHR*HRPRDA*TNSTRK 
JTRO  *  ATRUCK 
JFLO  *  BFLO 
RETURN 
END 

tENTRY  CIRSFA 
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